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Abstract. We present the results of ROSAT and XMM-Newton observations of the recurrent ultraluminous X-ray source 
(ULX) NGC 253 ULX1. This transient is one of the few ULXs that was detected during several outbursts. The luminosity 
reached 1.4xl0 39 erg s -1 and 0.5xl0 39 erg s -1 in the detections by ROSAT and XMM-Newton, respectively, indicating a black 
hole X-ray binary (BHXRB) with a mass of the compact object of > 1 1 M Q . In the ROSAT detection NGC 253 ULX1 showed 
significant variability, whereas the luminosity was constant in the detection from XMM-Newton. The XMM-Newton EPIC 
spectra are well-fit by a bremsstrahlung model (kT= 2.24 keV, Nh = 1.74xl0 20 cm -2 ), which can be used to describe a 
comptonized plasma. No counterpart was detected in the optical I, R, B, NUV and FUV bands to limits of 22.9, 24.2, 24.3, 22 
and 23 mag, respectively, pointing at a XRB with a low mass companion. 
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1. Introduction 

Ultra-luminous X-ray sources (ULXs) are extra-nuclear com- 
pact X-ray sources with luminosities considerably exceeding 
the Eddington luminos ity for stellar ma ss X-ray binaries of 
~ 2xl0 38 erg s" 1 (Makis hima et alfe OOO). However their lumi- 
nosities are still lower than that of active galactic nuclei (AGN). 

There are currently four preferred models to explain the lu- 
minosities of these objects. The first is that ULXs are inter- 
mediate mass black holes (IMBHs: M B h ~ 10 2 - 10 5 M G ). 
However IMBHs are at present not explainable with stellar 
evolutionary models. The alternatives are stellar-mass black 
hole X-ray binaries where either p hoton bubble ins tabilities al- 
low super-Eddington luminosities <Begelman2002), anisotrop- 
ically emitting X-ray binaries jKing et all200ll) . or that ULXs 
are micro-quasa rs that are observed do wn the beam of their rel- 
ativistic jet (e.g. lRevnolds et alJl997l) . It is therefore important 
to increase the sample of ULXs to find arguments that favour 
or exclude the above models. One attempt was the search for 
U LXs in 313 nearby ga laxies from ROSAT HRI observations 
bv lLiu & Bregmanl d2005l hereafter LB2005). A target of this 
search was the starburst galaxy NGC 253 where they found 21 
X-ray sources but only one of them matched their criteria for 
an ULX (NGC 253 ULX1). This source is located within, but 
close to the north-east boundary of the D25 ellipse of NGC 253. 

We here report on a more detailed analysis of 
NGC 253 ULX1 including ROSAT, XMM-Newton and 



* Based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA 



Chandra data, and specifically on the detection of a second 
outburst in one of the XMM-Newton observations. 



2. Search for the source in XMM-Newton, Chandra 
and ROSAT archives 

We searched the ROSAT, Chandra and XMM-Newton archive 
for observations of NGC 253. The results are listed in Table 
^ Except for two Chandra observations the position of 
NGC 253 ULX1 was always in the field of view (FOV). 
Besides the first detection in ROSAT observation 601 1 1 lh 
(LB2005), NGC 253 ULX1 was only visible in XMM-Newton 
observation 0110900101. These XMM-Newton and ROSAT 
HRI data are further discussed in Sect. 3 and 4, respectively. 

For the remaining observations we determined 3cr up- 
per limits for the count rate. From that we obtained up- 
per limit for fluxes and luminosities (cf Table 0. We used 
WebPIMMS (v3.6c) with the spectral model we got from the 
analysis of XMM-Newton observation 0110900101 to deter- 
mine energy conversion factors. The long term light curve of 
NGC 253 ULX1 is shown in Fig.|2] 

3. Detailed analysis of XMM-Newton observation 
0110900101 

NGC 253 ULX1 was detected for the second time on 2000 
December 14 with XMM-Newton. The position of the source 
was within the FOV of both of the MOS and the PN cam- 
eras. We used the latest version of the Science Analysis System 
(SAS v6.1) to process the obtained data from the PN (thin fil- 
ter) and the two MOS (medium filter) detectors. Periods of high 
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Table 1. Individual observations of NGC 253 ULX1. 



Date 


Instrument 


Observation ID 


Duration 


Flux" 


W 








(ks) 


/ -2 -1\ 

(erg cm s ) 


(yW erg s ') 


1991-12-08 


ROSAT 


600088h-0 


3.1 


< 1.1 x 10 -13 


< 0.08 


1991-12-25 


ROSAT 


600087p-0 


11.6 


< 2.2 x 10~ 14 


< 0.02 


1992-06-03 


ROSAT 


600087p-l 


11.2 


< 4.4 x 10~ 14 


< 0.03 


1992-06-05 


ROSAT 


600088h-l 


25.7 


< 4.4 x 10~ 14 


< 0.03 


1995-01-03 


ROSAT 


6007 14h 


11.0 


< 5.2 x 10~ 14 


< 0.04 


1 00£ AA 1 1 


n/OC A T 

KUbAl 


oUU / 14n-l 


in o 


< 2.U X tU 


< U.U2 


1997-12-20 


ROSAT 


601111h 


17.5 


1.8 x 10~ 12 


1.4 


1998-07-01 


ROSAT 


601113h 


2.0 


<2.8x 10 L3 


< 0.2 


1999-12-16 


Chandra 


969 




not in FOV 




1999-12-27 


Chandra 


790 




not in FOV 




2000-06-03 


XMM-Newton 


0125960101 


39.2 


<4.3 x 10~ 14 


<0.03 


2000-06-04 


XMM-Newton 


0125960201 


14.2 


< 1.8 x 10~ 14 


<0.01 


2000-08-16 


Chandra 


383 


2.16 


< 8.8 x 10 L5 


< 0.007 


2000-12-14 


XMM-Newton 


0110900101 


29.6 


6.3 x 10-' 3 


0.5 


2003-06-19 


XMM-Newton 


0152020101 


110.5 


< 3.6 x 10~ 15 


< 0.003 


2003-09-19 


Chandra 


3931 


83.6 


< 4.4 x 10~ 15 


< 0.003 



" 0.3-10 keV luminosity assuming a distance of 2.58 Mpc and a bremsstrahlung model (kT= 2.24 keV, N H = 1 .74x 10 20 cm 2 ) 



background were determined and excluded from further anal- 
ysis. The low background exposure times for PN and for the 
MOS instruments were 23.0 ks and 24.4 ks, respectively. 

We applied the source detection tasks eboxdetect and 
eraldetect only on the data from the PN detector, as the source 
was positioned far from the optical axis and close to the edge 
of the FOV on the MOS detectors. The obtained position was 
then corrected using optical reference coordinates from the 
USNO Bl catalogue JMonet et all2003l) of three AGN, identi- 
fied bv IVogler & Pietschl i 1 999t sources X4, X22, X58). The 
corrected position in J2000 coordinates is a — 00h48m20.11s, 
6 = -25° 10' 1074 with an error in position of 073. The derived 
position is well within the positional errors given by LB2005 
forNGC 253ULXl. 

A foreground star with a B magnitude of ~ 1 3 JMonet et alJ 
120031) is located close (15.5") to the obtained position (Fig.QJ- 
We can rule out that the actual detection of NGC 253 ULX1 in 
observation 0110900101 was caused by this star, as its proper 
moti on of -9.2 mas/yr i n RA- cos(DEC) and -3.6 mas/yrin DEC 
JZacharias et alJEo04l) is too small to match the detected posi- 
tion of NGC 253 ULX1 with that of the star within the period of 
observations. Additionally there was no detection of the source 
in other XMM-Newton observations using the same filter. 

We extracted energy spectra for NGC 253 ULX1 for all 
EPIC detectors. For the PN chip we included source counts 
from an elliptical region with major and minor axes of 27.6" 
and 12.3" respectively. The background region was a circular 
source-free region with a radius of 48" on the same CCD close 
to the source. For MOS the source extraction region was an 
ellipse with major (minor) axes of 28.95" (11.3") for MOS1 
and 31.35" (15.15") for MOS2, respectively. The background 
regions were circles with radii of 68" and 80" for MOS1 
and MOS2, respectively. After subtracting the background the 
spectra for each instrument were rebinned to a significance 
level of 3<r. 




0.5 1 2 5 

Energy (keV) 



Fig. 1. Comparison of the PN and MOS spectra of 
NGC 253 ULX1 with the best-fit bremsstrahlung model. In 
the lower panel the residuals (in units of <x) between data and 
model are shown. 

Table 2. Models for the source spectrum of NGC 253 ULX1 



Model" 


Nh" 
(10 20 cm" 2 ) 


Model parameter 


Xred 


po 


1.30 


T= 1.94 ± 0.05 


1.663 


bremss 
mekal 


1 74+0.02 
1.30 


kT= 2.24+°;^ keV 
kT=3.17±0.19keV 


0.961 
2.042 


diskbb 


1.30 


kT= 0.62 ± 0.04 keV 


1.671 


diskpn 


1.30 


kT= 0.69!°°* keV 
R=5.34tlHR s 


1.366 



8 po: power law, bremss: thermal bremsstrahlung, mekal: Mekal 
thermal plasma, diskbb: multiple blackbody disk, diskpn: accretion 
disk around a black hole 

b all models are modified by foreground absorption (XSPEC model 
wabs). 
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Julian Date - 2400000.5 

Fig. 2. Light curve of NGC 253 ULX1. Upper panel: solid 
symbols represent detections, open symbols 3<x upper limits of 
NGC 253 ULX1. Different instruments are represented by dif- 
ferent symbols: ROSAT (squares), XMM-Newton (diamonds) 
and Chandra (triangles). Lower panel: Single ROSAT HRI ex- 
posures with errorbars from observation 601 1 1 lh where the 
source was detected. The length of each observation is indi- 
cated by the x-errorbar. In contrast to the upper panel the lower 
panel plot is linear in luminosity. 



For the spectral analysis XSPEC 11.3.1 was used. The best- 
fit parameters from different models provided within XSPEC 
are listed in TablefJ] Using the PN and MOS spectra simultane- 
ously the source spectrum was best fitted with a bremsstrahlung 
model (Fig.[0. The fit of the multicolour disk blackbody model 
(diskpn) would also be acceptable. However, we favour the 
bremsstrahlung model since it is less complex and gives a bet- 
ter x 1 ,. Except for the bremsstrahlung model the foreground 
absorption (Nh) had to be fixed to the Gala ctic foreground ab- 
sorptionas alowerlimit(1.30xl0 20 crrT 2 , Dic kev & Lockmanl 
Il990l) . If the parameter was free to adjust it converged to un- 
reasonably low values. 

From the best fitting spectral model we calculated the 
sourc e flux and, assuming a distance of 2.58 Mpc JPuche et alJ 
119911) we derived an unabsorbed luminosity of 5. Ox 10 erg s 
in the 0.3-10.0 keV band. 

In order to study the temporal behaviour of the source a 
background corrected light curve was created using the tasks 
evselect and lccorr. The source count rate was constant at 
about 0.8 ct s within the errors of about 15% during observa- 
tion 01 10900101. 



4. Analysis of ROSAT observation 601 1 1 1 h 

The first detection of NGC 253 ULX1 was in ROSAT observa- 
tion 601 1 1 lh (LB2005). The observation (total exposure time 
17.5 ks) is spread over ten observing intervals, with different 
exposure and waiting time for the individual observations. The 
source was bright enough to determine luminosities for each of 
these observation intervals. 

We calculated count rates using the EXSAS source detec- 
tion task detect/sources. To reduce noise we only anal- 




Fig. 3. R-Band optical image taken with the Wide Field Imager 
on the MPG-ESO 2.2m Telescope. The source is located close 
to a ~ 13 mag star. The R magnitude is 24.2. 




Fig. 4. Near UV image taken with GALEX. The NUV magni- 
tude is 22. 



ysed HRI channel 2-15. We used WebPIMMS (v3.6c) and the 
spectral model retrieved from the XMM-Newton observation 
(bremsstrahlung, kT= 2.24 keV, N H = 1.74xl0 20 cm" 2 ) to de- 
termine energy conversion factors to obtain the corresponding 
fluxes and luminosities (see lower panel of Fig.©. The lumi- 
nosity averaged over the whole observation (1.43xl0 39 erg s _1 ) 
is indicated by the dashed line. 

During the observation the source showed significant vari- 
ability by at least a factor of 2. 



5. Discussion 

We detected the recurrence of NGC 253 ULX1 in the XMM- 
Newton observation from 2000 December 14. This was the first 
detection after the outburst in 1997, reported from ROSAT HRI 
observations by LB2005. In all other observations of NGC 253 
the luminosity of the source was below the detection limit. This 
implies brightness variability by at least a factor of 500. Its 
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fastest change in luminosity (Lmax/L^jj) exceeds a factor of 
71 in 120 days. 

The improved position (errors of 0V24 compared to 4"- 
10") of the source determined in Sect. 3 allowed us to search 
for optical counterparts. We checked images taken with the 
Wide Field Imager (WFI) on the MPG-ESO 2.2m Telescope 
at La Silla in the R- (Fig. [3j, I- and B-band (limiting magni- 
tudes 24.2, 22.9 and 24.3, respectively) and images taken with 
the Galaxy Evolution Explorer (GALEX, a space telescope 
from NASA observing in the ultraviolet) in the NUV (Fig. |4j 
and FUV (limiting magnitudes 22 and 23, respectively), but no 
counterpart could be detected. 

With the data discussed in Sect. 3 and 4 we can exclude 
that NGC 253 ULX1 is either a foreground object or a back- 
ground AGN based on three arguments: (i) We estimated the 
log(/x//opt) = log/x + {my 12.5) + 5.37 > 3.2 using the 
flux of the ROSAT detection and a lower limit for my of 
24.2 (averaging the limiting magnitu des of the R- and the B- 
WFI images, see above). Following iMaccacaro et al] fl988) 
this value exceeds that expected for galactic sources (-4.6 to 
-0.6) as well as AGNs (-1.2 to +1.2). (ii) The variability of 
NGC 253 ULX1 is by a large factor higher than the typical 
value observed for AGNs (~ 10 - 60). (iii) NGC 253 ULX1 
shows a bremsstrahlung spectrum, whereas spectra of AGNs 
above 2keV are typically fitted by a power law. The recur- 
rent outbursts also exclude that the source is the luminous 
remnant of a recent supernova, li ke e.g. SN1993J in M81 
(Zimmermann & Aschenbach 2003). 

The X-ray spectrum may indicate that NGC 253 ULX1 
is a low mass X-ray binary (LMXB). The X-ray emission in 
these objects is created in the optically thin boundary layer 
between the disk and the neutron star and compton ization 
may dominate the spectral emission llWhiteetalj fl988i. lead- 
ing to a spectrum that can be fitted with a bremsstrahlung 
model. However the ROSAT HRI peak luminosity of 1.43 
xlO 39 erg s _1 is very high for typical LMXBs. Other systems 
that show bremsstrahlu ng spectra are black hole XRBs, e.g. 
Cv gX-1 JSunvaev & Truemperlll979h . LMC X-3 and X1755- 
33 ffWhite et alJll988l) . These systems may contain a high or 
low mass companion. 

An additional argument for a low mass companion comes 
from the lack of an optical counterpart (see above). High mass 
X-ray binaries (HMXBs) should be detectable at about 22 
to 24 mag, extrapol ating V magnitudes from HMXBs in the 
Magellanic Clouds fciu et alJl2000h . We would have detected 
an object of this brightness in the WFI data. 

The luminosity of a compact object radiating at the 
Eddington limit is given as L El u = 1.5 X 1O 38 (M/M ) erg s _1 , 
when electron scattering dominates the opacity. Luminosities 
higher than 2xl0 38 erg s _1 (corresponding to a 1.4M object, 
commonly assumed as the maximum mass of a neutron star) 
suggests that the compact object is a black hole. According to 
NGC 253 ULXl's maximum luminosity of 1.43xl0 39 erg s" 1 
the lower limit for the mass of the black hole is 11 M . 
Therefore NGC 253 ULX1 is not required to be an IMBH. 

Another argument against an IMBH is the temperature of 
NGC 253 ULX1. Mill er et alJ (Hp04a) compared intermediate 
mass black hole candidate ULXs and stellar mass black holes 



with respect to luminosity and temperature. If we assume the 
multicolour disk blackbody model then NGC 253 ULXl's po- 
sition in the lumin osity-disk temperature diagram (Fig. 2 in 
iMiller et alJl2004al) indicates that NGC 253 ULX1 is not an 
IMBH, but a stellar mass black hole. 

Recently another object was found that, like 
NGC 253 ULX1, showed also a bremsstrahlung spec- 
trum: X-44 in the Antennae Galaxies (NGC 4038/4039) 
(Mil ler etal.l2004bl) . The temperature of X-44 is 3.7 + 0.5 keV 
and its luminosity is 1 -0*i - | XlO 40 ergs" 1 . This temperature 
is about a factor of 1.5 higher than in NGC 253 ULX1, and 
the luminosity exceeds the luminosity of NGC 253 ULX1 by a 
factor of 15 compared to the outburst in 1997. 

Another inter esting ULX to com pare NGC 253 ULX1 with 
is M101 ULX-1 jKong et alJl2005l) . It was the first ULX that 
like NGC 253 ULX1 has been observed during more than one 
ultra-luminous outburst. Like many other ULXs the spectrum 
of M101 ULX-1 is best described with an absorbed black- 
body model, but the temperature of ~ 50 - 160 eV is rather 
low. M101 ULX-1 has a peak luminosity of about 10 41 erg 
s _1 (0.3 - 7keV), and the hardness of its spectrum changed 
between different observations. We do not know whether the 
spectrum of NGC 253 ULX1 changed in the two observations, 
as the XMM-Newton data provided the very first spectrum of 
the source. During 12 years of observations NGC 253 ULX1 
showed two outbursts with an interval of three years. In 
M101 ULX-1 the two outbursts are only seperated by half a 
year. On shorter time scales NGC 253 ULX1 showed only one 
drop in luminosity by a factor of ~ 2 during the ROSAT obser- 
vation, and in the XMM-Newton observation (exposure time 
8.2 h) no variability could be detected. M101 ULX-1 on the 
other hand does show short-time-scale variability. Its luminos- 
ity changed by a factor of > 10 on a time scale of hours. The 
lack of short time variability of NGC 253 ULX1 argues against 
the relativ istic beam ing model, since this would require a very 
stable jet jRevnolds et alJl997l) . 
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